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The coexistence of the universal transverse force and the non-universal friction force on a topo- 
logical singularity in a superfiuid is shown here. Based on the BCS type microscopic theory, we 
explicitly evaluate the quasiparticle contribution to the friction coefficient in a clean fermionic su- 
perfiuid, showing a new feature of logarithmic divergence. 
PACS numbers: 67.40.Vs; 47.37.+8 



I. INTRODUCTION 

When a vortex in a superfiuid moves, in addition to the hydrodynamic Magnus force which 
is proportional to the superfiuid density, the coupling of the vortex to quasiparticles, phonons and 
impurities may cause extra forces in both transverse and longitudinal directions of the vortex motion. 
Those possible non-vanishing forces have been the underlying assumption for phenomenological 
models used in analyzing experimentsEI, though analyses based on topological methods such as the 
Berry phase always give the universal transverse force.0 

While the total transverse force on a moving vortex needs further investigations, the friction force 
is poorly studied. The usual quoted friction coefficient formula derived from a microscopic Hamil-i 
tonian is obtained through a relaxation time approximation in force-force correlation functions^ 
However, this type of procedure is incorrectu, for the reason known since the 60's in the context 
of obtaining the friction of electrons using force-force correlation functions. The lack of an explicit 
and correct microscopic derivation of the friction has two consequences. First, though the friction of 
vortex motion is an important quantity in many experiments, it has not been received due attention. 
Second, the incomplete understanding of the friction has generated doubts about the exact result 
on the transverse force. 

Here we describe a self-contained theory of both the friction and the transverse forces by a 
rigorous and elementary method. The method we choose follows that of Ref.[5]. In order to obtain 
the friction, we use a limiting process similar to the one in transport theory and in non-equilibrium 
statistical mechanics, which has not been explicitly discussed in Ref. [5]. 

In the following, after giving a general expression for the friction in Section 2, as an example, we 
evaluate the friction for the case of a clean fermionic superfiuid in Section 3. The resulting friction 
coefficient is new. It comes from the off-diagonal potential scattering of the extended quasiparticles, 
and is stronger than the ohmic damping by a logarithmic diverging factor. 



II. TRANSVERSE AND FRICTION FORCES 



We consider an isolated vortex in the superfiuid whose position ro is specified by a pinning poten- 
tial. The system is otherwise homogeneous and infinite. There is no externally applied supercurrent 
or normal current. The vortex is allowed to move slowly. Hence the system Hamiltonian H contains 
a slowly varying parameter ro(t). The many-body wavefunction of the superfiuid can be 

expanded in terms of the instantaneous eigenvalues _Eo,(ro) and eigenstates \ip a (ro)), for which we 
choose phases such that {ip a \ipa) = 0. Because of our assumption of a homogeneous superfiuid, 
E a (ro) is independent of both ro and time t. With those considerations, the many-body wavefunc- 
tion \^ a [t)) can be expressed by 

|*«(*)) = e- iE " t/n \^ a {r )) + a a ,(t)e- lE °' t/h \iP a ,(r )} , (1) 

a' T^a 

where, to the first order in velocity, a Q (t) = 1, and 

a a ,{t) = - [ dt'<^|^)e l(BQ '" Ba)t ' /R ■ (2) 
Jo 
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This gives the expectation value of the force on the vortex as 



F = -^/ Q (* Q |V ^|*a> 

a 

= -^f a {^ a \V H\^ a ) + J2 /c«(^(ro(t))|Votf|i/vMt))} x 

t 

dt\i; a ,(r (t'))\MMt'))}e i(E "'- Ea)(t '- t)/h + c.c. , 

where f a is the occupation probability of the state a. For a vortex moving with a small and uniform 
velocity vy = ro(t), 

\ip a (ro)) = vy • \Voipa(ro)) ■ 

Here Vo denotes the partial derivative with respect to the position ro of the pinning potential, the 
vortex position. The first term in the right hand side of Eq.(3) is independent of vv and will be 
ignored. In fact, it is zero because of the translational invariance regarding the vortex position. 
The integration in time can be carried out directly. In addition, we also need to use the following 
relations, 

GMVoiWa'} = (E a > - E a ){il) a \Voi> a ,) = (E a - S^XVoValVV') , 

which are obtained by taking gradient Vo with respect to H\ip a i) — E a i\ip a i) and (lp a \H = E a {ip a \, 
then multiplying from left or right by {ip a \ or |Vv) respectively. 

To obtain the long time behavior of Eq.(3) requires a limiting procedure. There are two ways 
of taking limiting sequences. A possible sequence is to take the low frequency limit before the 
thermodynamic limit. As we will find out, the transverse force is independent of the limiting process. 
Therefore such a limiting sequence is correct as long as only the transverse force is concerned. In 
such a calculation, because the energy levels have been treated as discrete ones, there is no friction. 
In order to obtain friction, we have to take the thermodynamic limit before the low frequency 
limit. For this purpose, we use the Laplace averagecl, limt^oo F(i) = lim e _> + F(e), where F(e) = 
e J °°di F(t)e~ £ ', and F(e = 0) = Fj_+F||. We will also use the identity l/(e+ix) = ir5(x)-iP(l/x). 
For the transverse force Fj_, we have 

lim e / ctoT £ ' (1 - cos(£ a - E a ,)t/h) = lim (1 - irS(E a - E a ,)e) = 1 , 

because the summation over states in Eq.(3) is well behaved regardless whether E a is discrete or 
continuous. Hence 

Fi=ffi)^ f a {((V>a|VoVv) x {V ip a i\ip a }) ■ z} w X z . 

a' -£cl 

This is precisely what has been obtained in Ref.[5], where further calculations lead to Fx = 
—hLp s vv x z, independent of details of the system. Here p s is the superfluid number density, 
L the length of the vortex. The-, universal nature of the transverse force has been confirmed ki 
experiments of the vibrating wired and the vortex precessiona, and even in dirty superconductors!! 
The longitudinal force, friction, is given by Fy = — r)Vv with 

= \ E h E a 7- f E a 5{Ea ~ ^')K^|VoH|^)| 2 . 

a' -£cl 

The friction coefficient 77 is determined by low energy excitations such as phonons, extended quasi- 
particles, and localized quasiparticles when their discrete energy spectrum is smeared out by im- 
purities. This expression is identical to the result in Ref. [10] for the case of ohmic damping in the 
zero frequency limit. Eq.(7) will not pick up any superohmic contributions, and will give infinity for 
any subohmic contributions. We point out that with the aid of Eq.(5) the transverse force can be 
expressed only in terms of the wavefunction or the density matrix without explicit referring to the 
Hamiltonian or its eigenvalues, as shown by Eq.(6), on the other hand, the longitudinal force, the 
friction, cannot, as shown by Eq.(7). The explicit dependence on the Hamiltonian or its eigenvalues 
in Eq.(7) is the source of the sensitivity of friction to details of the system. 



III. QUASIPARTICLE CONTRIBUTION 



For the superfluid 4 He, there is no full microscopic theory yelE. We will not attempt to evaluate 
Eq.(7) for this superfluid here because of the sensitivity of the friction to details. The situation can 
be very different in the case of the recent Bose-Einstein condensed systems, where a well defined 
microscopic theory is supposed to be known. Instead, as an example to illustrate the directness 
and usefulness, we evaluate Eq.(7) for the case of a homogeneous fermionic superfluid using BCS 
theory with s-wave pairing. At finite temperatures the extended states above the Fermi level, the 
quasiparticles, are partially occupied. The vortex motion causes transitions between these states 
and gives rise to friction. The transitions between different single quasiparticle levels (ij) a \VoH\il) a i) 
are considered here since they dominate the low energy process. The quasiparticles are described by 
the eigenstates, u a and v a , of the Bogoliubov-de Gennes equation. Their behavior in the presence 
of a vortex has been well studied in Ref.[ll]. We may take 

^ = ( ZU ) = ^e^e^ + ^ 2 f(r) , (8) 

with r measured from the vortex position, and 9 the azimuthal angle around the vortex. In order 
to obtain an analytical form for the transition element, we use a WKB solution for ,f(r), 



[1 ± y /E* - \Hr)\ 2 /Ef/ 2 
2/r7 V [1 =F V E2 -|A(r)| 2 /,E] 1/2 



cxp 




dr' (k 2 p r -^±±^^/Ei-\A(r>W) \ + c.c. . (9) 
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Here k 2 = k 2 t — k 2 , R is the radial size of the system. This WKB solution is valid when r is outside 
the classical turning point r t = \fj,\/k p . Here r t is the impact parameter. A WKB solution also exists 
inside the turning point. However, because it approaches zero as (rfc p )' M '/|^|!, the contribution to 
the transition elements from this region is small, and will be set to zero. The transition elements 
are then given by 



|<vgVo#|W>| 



dx(u* a ,(x)VoAv a (x) + Va'(x)S7oA*u a (x)) 



, > £ofc p 

Here A^ is the value of | A(r)| far away from the vortex core. Physically it means that if the classical 
quasiparticle trajectory is far away from the vortex core, it will not contribute to the friction. The 
summation over states in Eq.(7) is replaced by 



E 



Ol'^Ol fix ,t^2, k zl, k z2 



dEld E 2 - E ^ E * I m ^ R 



y/E* - A 2 , y/E£ - Alo \h z k f 



after considering the density of states. 

Substituting Eq.(lO) into Eq.(7), using the quasiparticle distribution function f a — l/(e l3Ea + 1), 
the coefficient of friction is given by 

Lm%AtP f, p E 2 1 . 

V Stt^ 3 J Aq E* - AL £2 cosh 2 (PE/2) ' ( ' 

The integral in Eq.(ll) diverges logarithmically. It implies that the spectral function corresponding 
to the vortex-quasipaptieie coupling is not strictly ohmic but has an extra frequency factor propor- 
tional to ln(Aoo/?iu;)L3il3. When Tllu is not very small comparing to Aoo, which may be realized 
when close to T c , we can ignore the logarithmic divergence in Eq.(ll) by using the density of states 
for normal electrons to obtain a finite friction, i.e. replacing E 2 /{E 2 — A 2 ^) with 1 in Eq.(ll). Close 
to T c , the friction approaches zero the same way as A^, which is proportional to the superfluid 
density p s . When — ln(7iw/Aoo) is large, we need to use a more accurate expression of vortex friction 
obtained in Ref. [10]. Straightforward evaluation shows that in such a case 



Lm 2 CoAL/3_ 

"cosh 2 (/3Aoo/2) 



V = — 2 ,L M MAoo/ft^) . (12) 



3 



Here u> c is the low frequency cut-off. ft is determined by the size of the system for a single vortex, 
and by the inter-vortex distance for a vortex array. 

It should be emphasized that the logarithmic divergence comes from the interplay between the 
divergence in the density of states and the off-diagonal potential scattering. We can consider a 
situation in which we physically create a pinning center to trap the vortex and guide its motion. 
In such a case the vortex has a diagonal potential. If the scattering is dominated by the diagonal 
potential, e.g., by the pinning potential, an additional factor coming from \u a \ 2 — \v a \ 2 will remove 
this logarithmic divergence. This again shows the sensitivity of the friction to system details. 

The above results apply to the dynamics of the axisymmetric vortex in 3 He B phase, where the 
Bogoliubov-de Gennes equation is essentially the same as the one in a s-wave superconductor. 



IV. CONCLUSION 



The transverse force formula confirms the Berry phase results obtained with a trial many-body 
wavefunction. New results on the friction have been obtained from quasiparticle contributions, 
which can be further tested experimentally. 
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